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Abstract 

The bifurcation structure of the Langford equation is studied numerically 
in detail. Periodic, doubly-periodic, and chaotic solutions and the routes to 
chaos via coexistence of double periodicity and period-doubling bifurcations 
are found by the Poincarc plot of successive maxima of the first mode Xi. 
Frequency-locked periodic solutions corresponding to the Farey sequence F„ 
are examined up to n = 14. Period-doubling bifurcations appears on some of 
the periodic solutions and the similarity of bifurcation structures between the 
sine-circle map and the Langford equation is shown. A method to construct 
the Poincare section for triple periodicity is proposed. 
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1 Introduction 



A quasi-periodicity route to turbulence due to Ruelle and Takens [7] is that transition could 
occur via three successive Hopf bifurcations, leading from a fixed point to a limit cycle, then 
to a torus and finally to a 3-torus. On the other hand, double periodicity is considered to be 
well modeled by the one dimensional sine-circle map, which indicates that frequency-locking 
and periodic solutions appear as the nonlinear parameter increases. Therefore, it is a natural 
question what actually happens on attractors of a simple set of ordinary differential equations 
(ODEs) having double periodicity, like the Langford equation [4-6]. Similar numerical studies 
have been done for the five dimensional ODEs modeling magnetoconvcction [1, 2] and the six 
dimensional ODEs of the Gledzer shell model of turbulence [8] . In the latter case the parameter 
to be changed is the viscosity, or equivalently, the Reynolds number. Both studies show a 
bifurcation structure very similar to that of the sine-circle map, although triple periodicity is 
stated based on the numerically obtained Lyapunov exponents in [2]. The presented results 
correspond to a detailed study that extends Langford [4-6]. 

In Section 2, a method to construct the three dimensional Poincare section is proposed. For 
triple periodicity expected by Ruelle and Takens [7], points lie on a surface. Frequency- locked 
double periodicity indicates points on a closed curve embedded on the surface. Section 3 gives the 
Langford equation with an explanation of the evolution of the energy and selected parameters. In 
Section 4, numerical results of bifurcation structures of the equation are shown. It is confirmed 
that, instead of triple periodicity, we have double periodicity, frequency-locking and period- 
doubling bifurcations. The structure is very complicated; frequency-locking corresponding to 
the Farey sequence with the index up to 14 is confirmed, and plural sequences of the period- 
doubling bifurcations are observed. Summary and further possibilities to explore complexity 
of dynamical systems modeling fluid dynamics and other high dimensional systems with the 
presented approach are described in Conclusions. 

2 Triple periodicity 

As a typical example of triple periodicity, we consider the following function which is a sum of 
three sinusoidal functions of t; 
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x{t;a,uj) = sint + sinnt + asinujt. (1) 

A standard method to construct the Poincare map used in [1,8] is to seek for successive local 
maxima x„ of x{t). Then we let = Xn + ixn+i and 9n = (27r)~^Arg(z„) (mod 1). For double 
periodicity, 6n obeys the generalized sine-circle map [8], which gives a curve on the {6n,0n+i) 
plane. Similarly, we can anticipate a surface in the {On, On+i, space. To check this analogy, 

we show the Poincare plot of Eq. ([1]) in the time interval < t < 5000 for double periodicity 
a = in Figure 1, triple periodicity a = 0.5, u; = \/2 in Figure 2, and frequency-locked double 
periodicity a = 0.5, u; = 1.4 in Figure 3. Points lie on a curve in Figure 1 (a) and on a torus in 
Figure 2 (a), whose three dimensional structure can be viewed by a new 3D graphics function in 
Mathematica 6. The latter indicates a possibility that 9n+2 can be expressed by a function of On 
and On+i- This Poincare plot will give a method to identify triple periodicity among complicated 
time series of general numerical or observed data. 

We have 5:7 frequency-locking in the case of Figure 3. Correspondingly, points are again 
on a curve, which is more complicated than that in Figure 1. Actual data may contain higher 
harmonics in Eq. ([T]) but the qualitative behavior can be expected to be similar. 

3 The Langford Equation 

The Langford equation is the set of three ordinary differential equations for Xi{t), i = 1,2,3, 
given as follows: 

{X3 - b)xi - CX2, (2) 

CXi + (X3 - b)x2, (3) 
xl 

d + ax3 - — 

-{xl + xl){l + fx3) + ex3xl (4) 
The temporal evolution of the energy defined by = (1/2) X]?=i is given by 

^ =dx3 + axl-^ -{b + fxl){xl + xl) + ex^xl (5) 



1^= = 

dX2 

-dr= = 

dX3 
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The set of the parameters fixed in this paper is borrowed from [6] as 



(a, b, c, d, f) = (1, 0.7, 3.5, 0.6, 0.25), (6) 

and e is the changing parameter. Because 6 > and / > 0, the right hand side of ([5]) becomes 
negative if E is sufficiently large for e = 0, leading to the proof that the solution is finite. If 
e 7^ 0, the finiteness of the solution is unknown but supported by the numerical solution for 
small e. 

The contraction rate of the volume element in the phase space is given by 

dxi dxo dxs 9 , „/ 9 9x /_s 

IT^ + ir^ + ir^ = -4 + ^^3-^b + a- f{xl + xl) + ex?. (7 
oxi 0x2 0x3 

It depends on the position in the phase space. If e = 0, it becomes negative for sufficiently large 
values of E since / is positive and the quadratic terms become dominant. 

4 Bifurcation Structure 

Before we show numerical results of the Langford equation, it is useful to review the bifurcation 
structure of the sine-circle map: 

9n+i = f{9n) = 9n + n + ^sm{27T9n), (mod 1). (8) 

zvr 

In order to compare the diagram of ([8]) with that of the Langford equation, the parameter fl 
and K are parametrized as 

= 0.5(1 + a), K = ia. (9) 

The initial condition is ^1 = 0, < a < 1, Aa = 0.001, the total iteration is 400, and the 
last 200 steps are plotted in Figure 4. Although double periodicity and frequency-locking are 
illustrated by the sine-circle map in many literatures, the period-doubling bifurcation is also 
observed at a ~ 0.6 in Figure 4, which can make it easy to understand the similarity of the 
bifurcation structures between the sine-circle map and the Langford equation. The coexistence 
of two scenarios of routes to turbulence, quasi periodicity and period-doubling bifurcation, may 
be common in many of dynamical systems. 



In order to study bifurcation structures of the Langford equation, sets of xi at its local 
maximum after transient states are plotted with various values of the parameter e. Eqs. ([2][5]) are 
solved by NDsolve command in Mathematica 5.2, which makes 100 to 500 different computations 
possible in a single program. The typical final time is t/ = 2000 ~ 4000 and the last period 
tf — td < t < tf with td = 200 ~ 600 is picked up to identify attractors. In order to find 
the local maximum, the time period is divided into intervals with the width At = 0.1, and 
for the interval which can include the local maximum, the FindMaximum command is invoked. 
The typical initial condition for the numerical computation is {xi,X2,X3) = (0.01,0,0). We 
also choose the initial condition as the final state just before the parameter varied in order to 
examine the hysteresis. In some cases, multiple stable states are observed. About 40 different 
runs are performed with various regions of e and suitable numerical parameters. 

Figure 5 shows the bifurcation diagram for the parameter range < e < 0.2. The step of 
the parameter e is Ae = 0.0001. The region e < 0.03 includes doubly periodic solutions. Then 
3:4 frequency-locking yields the periodic solution. The Feigenbaum period-doubling bifurcation 
occurs at e ~ 0.057. There are many periodic windows as well as chaotic solutions after the 
period-doubling bifurcation. Its behavior is very similar to that of the one dimensional logistic 
map. 

The range 0.02 < e < 0.03 is enlarged in Figure 6, in order to see the most dominant 
frequency-locked periodic solutions. The solution is doubly periodic if e < 0.02028, even if e 
is negative. As e increases, a periodic solution with 17:23 resonance appears at e ~ 0.02028. 
Periodic solutions shown by m2 points are specified by the resonance condition rrii : m2 drawn 
at just below the upper frame of Figure 6. The region including double periodicity ends with 
the emergence of 3:4 resonance at e ~ 0.02935. Doubly periodic solutions draw points whose 
numbers are controlled by the time range for plotting and typically much larger than the periodic 
solutions. 

The Farey sequence for integers n > is the set of irreducible rational numbers a/b where 
< a < b < n and the greatest common divisor of a and 5 is 1. The Farey sequence appears in 
frequency-locking of the sine-circle map. The resonance conditions between 17:23 and 3:4 are 
shown up to the Farey index n = 14 in Table 1. 

For the 23:31 resonant periodic solution, an incomplete period-doubling bifurcation is ob- 
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Tabic 1: Farey sequence up to n = 14. 




served in Figure 6. Figure 7 shows the diagram enlarged for 0.021 < e < 0.0235. Figure 8 shows 
the diagram for 0.0243 < e < 0.0253. The period-doubling bifurcations observed on both of the 
two parameter regions of 43:58 resonance. 

The periodicity can be also judged by numerical computation of the square of the minimum 
distance 
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D = mmY,{xi{t)-xm?, (10) 
i=l 

indicating the accuracy of recurrence. In Figure 9 for 0.0272 < e < 0.0283, the small value of 
D indicates frequency-locked periodicity including resonances corresponding to n = 15 (35:47) 
and n = 16 (38:51 and 93:125) of F„. 

The width of resonance decreases as the index n of F„ increases, that implies scaling laws 
as noted by [1]. Summary of observed stable periodic windows is given in Table 2. Figure 10 
shows the parameter e versus the Farey index n for the periodic windows. Many of resonances 
mi : 777,2 corresponding to Farcy sequence Fn are confirmed, although some of them are not 
observed because of the possible lack of stability. 
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5 Conclusions 



In conclusion, bifurcation structures similar to [1,2,8] is observed in the Langford equation; the 
coexistence of the double periodicity, frequency-locking and period-doubling bifurcations. The 
Langford equation would be one of the most illustrative ODE systems having double periodicity 
since it consists of only three variables. Each periodic attractor represents a corresponding 
knot in three dimensional space, as discussed in [3] for the Lorenz system and in [2] for the 
magnetoconvection system. 

Recent progress of cost performance in computer hardware would make it possible to judge 
whether similar double periodicity and frequency-locking appear in the high dimensional Navier- 
Stokes turbulence and other chaotic systems. The method to judge quasi periodicity would be 
basically the same as that stated in the presented work. 

This study was motivated by comments by Professor Bekki at the meeting of Japan Physical 
Society several years ago. The author is grateful to Professor Yamagata for support of research 
and to Professor Langford for sending the author his reprints. 
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Figure 1: Poincare plot of ([T]) for the double periodicity a = in (a) the (a;„, Xn+i, Xn+2) space, 
on (b) the {9n,9n+i)) plane, and in (c) the (0^, ^n+i, space. Points A, B and A', B' in 

Figures 1(a) and 1(b), and C, C in Figure 1(c), which are separated due to modulus 27r in the 
argument, are the same points. In order to see which surface the points lie on, the rectangles 
including the points are added in Figure 1(c). 
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Figure 2: Poincare plot of ([T|) for the triple periodicity a = 0.5, tj = \/2 in (a) the 
{xn,Xn+i,Xn+2) space, Oil (b) the {6n,6n+i) plane, and in (c) the {0n,9n+i,0n+2) space. The 
points are on the surface in Figure 2(a) and 2(c). As a in Eq. ([T|) increases, the width of the 
surface becomes large. 
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Figure 3: Poincare plot of ([T]) for the frequency-locked double periodicity a = 0.5, w = 1.4 in 

(a) the {xn-,Xn+i,Xn+2) spacc and (b) the 6'„+i, 6'„+2) space. 
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Figure 4: Bifurcation diagram of the sine-circle map ([8]). The parameters O and K are related 
as ([9]). 40000 points are randomly selected from the data for a clear view. 
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Figure 5: Bifurcation diagram of the Langford equation for < e < 0.2. 40000 points are 
randomly selected from the data for a clear view. In the following figures, bifurcation diagrams 
are for the Langford equation. 
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Figure 6: Enlarged bifurcation diagram for 0.02 < e < 0.03. 
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Figure 7: Bifurcation diagram enlarged again for 0.021 < e < 0.0235. 
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Figure 8: Enlarged bifurcation diagram for 0.0243 < e < 0.0253. Two parameter regions for 
43:58 and 83:112 resonances are observed. The period-doubling bifurcation follows from 43:58 
periodic solutions. 
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Table 2: Observed periodic windows, mi : m2 denotes the resonance condition, n the Farey 
index, and (e„) the lower (upper) limit of e, respectively. NF denotes that the resonant 
periodic solutions are not found. 



mi : m2 


n 






17:23 


9 


0.02028 


0.02112 


88:119 


14 


0.021625 


0.021640 


71:96 


13 


0.021794 


0.021817 


125:169 


14 


0.021924 


0.021926 


54:73 


12 


0.022056 


0.022113 


145:196 


14 


0.022217 


0.022219 


91:123 


13 


0.022287 


0.022299 


128:173 


14 


0.022373 


0.022377 


37:50 


11 


0.022513 


0.02269 


131:177 


14 


0.022813 


0.022822 


94:127 


13 


0.022873 


0.022898 
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0.022937 


0.022947 
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0.024605 
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0.0247 


43:58 
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0.02179 
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109:147 


14 


NF 
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